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Abstract — A new construction of a directional continuous 
wavelet analysis on the sphere is derived herein. We adopt 
the harmonic scaling idea for the spherical dilation operator 
recently proposed by Sanz et al. [1] but extend the analysis to a 
more general directional framework. Directional wavelets are a 
powerful extension that allow one to also probe oriented structure 
in the analysed function. Our spherical wavelet methodology has 
the advantage that all functions and operators are defined directly 
on the sphere. The construction of wavelets in our framework is 
demonstrated with an example. 

Index Terms — Wavelet transforms, spheres, convolution. 
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I. Introduction 

THE extension of wavelet analysis to a spherical manifold 
is of important theoretical and practical interest. Wavelets 
are a powerful signal analysis tool due to the spatial and scale 
localisation encoded in the analysis. The usefulness of such 
an analysis has previously been demonstrated on a large range 
of physical applications. However, many of these applications 
are restricted to Euclidean space R n , where the dimension 
of the space n is often one or two. Nevertheless, data are 
often defined on other manifolds, such as the 2-sphere (S* 2 ; 
hereafter we refer to the 2-sphere simply as the sphere). For 
example, applications where data are defined on the sphere 
are found in astrophysics (e.g. [2], [3]), planetary science 
(e.g. [4]-[6]), geophysics (e.g. [7]-[9]), computer vision (e.g. 
[10]) and quantum chemistry (e.g. [11], [12]). A wavelet 
transform on the sphere is required to realise the potential 
benefits of a wavelet analysis in these settings. In this work 
we construct a new directional continuous spherical wavelet 
' transform (CSWT). 

A number of attempts have been made previously to 
extend wavelets to the sphere. Discrete second generation 
wavelets on the sphere that are based on a multi-resolution 
analysis have been developed [13], [14]. Haar wavelets on 
the sphere for particular tessellation schemes have also been 
developed [15], [16]. However, these discrete constructions 
may not lead to a stable basis. Furthermore, the symmetry 
of the sphere is clearly lost when a tessellation scheme is 
adopted. Other authors have focused on the continuous wavelet 
transform on the sphere. Continuous spherical wavelets have 
been constructed directly from spherical harmonic functions or 
Legendre polynomials by [17], [18], however these solutions 
suffer from poor localisation and require an abstract dilation 
parameter that imposes a number of ad hoc assumptions. A 
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satisfactory extension of the continuous wavelet transform to 
the sphere is defined by [19], however this construction also 
requires an abstract dilation parameter and corresponding ad 
hoc assumptions. Others adopt a tangent bundle viewpoint 
[20], [21], thereby avoiding the necessity to define a dilation 
operator on the sphere. Many of these previous works also 
consider only azimuthally symmetric wavelets and do not 
make the extension to a directional analysis. A consistent 
and satisfactory framework for a directional spherical wavelet 
transform has been constructed and developed by [22]-[28]. 
This transform was originally derived entirely from group 
theoretic principles [22], however the formalism has since been 
reintroduced independently of the original group theoretic 
formalism in an equivalent, practical and consistent approach 
[28]. A stereographic transform is adopted in these works 
to relate spherical wavelets and operators to their planar 
counterparts. For a more detailed review of the attempts made 
to construct a wavelet transform on the sphere see [22], [23], 
[29]. 

Recently, an alternative wavelet construction on the sphere 
has been derived by [1]. In this setting the spherical di- 
lation operator is defined by a scaling operation in spher- 
ical harmonic space. This construction has the advantage 
that all functions and operations are defined directly on the 
sphere, rather than through projections to and from the plane. 
However, the derivation presented by [1] is not fully general 
and is restricted to azimuthally symmetric wavelets. We adopt 
the harmonic scaling idea for the dilation operator proposed 
by [1] but extend the spherical wavelet analysis to incorporate 
directional wavelets. Directional wavelets are a powerful ex- 
tension that allow one also to probe oriented structure in the 
analysed function (see e.g. [25] for a discussion of directional 
wavelets). 

The remainder of this paper is organised as follows. We 
begin by outlining some mathematical preliminaries in sec- 
tion II. In section III and section IV we derive the spherical 
wavelet analysis and synthesis respectively. In section V we 
show some properties of our spherical wavelet construction. 
We demonstrate in section VI the construction of spherical 
wavelets in our formalism. Finally, concluding remarks are 
made in section VII. 

II. Mathematical preliminaries 

It is necessary first to outline some mathematical prelimi- 
naries before constructing the wavelet transform on the sphere. 
By making all assumptions and definitions explicit we hope to 
avoid any confusion over the conventions adopted. The reader 
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familiar with harmonic analysis on the sphere may skip this 
section and refer back to it as required. 

We consider the space of square integrable 
functions L 2 (S 2 , dfl) on the unit 2-sphere S 2 , where 
dQ(co) = sintfdtfd^ is the usual rotation invariant measure 
on the sphere and uj = (#, G S 2 denotes spherical 
coordinates with colatitude 6 and longitude <p. A square 
integrable function on the sphere / G L 2 (S 2 , dfl) may be 
represented by the spherical harmonic expansion 

oo i 
£=0 m=-£ 

where the spherical harmonic coefficients are given by the 
usual projection on to the spherical harmonic basis functions 



The * denotes complex conjugation. We adopt the Condon- 
Shortley phase convention where the normalised spherical 
harmonics are defined by 



w«)-<-ir^^jr<~«.~, 

where PJ n {x) are the associated Legendre functions. Using 
this normalisation the orthogonality of the spherical harmonic 
functions is given by 



i^ m (a;)Y? m ,(a;) = fe/<Wi' , (1) 



where Sij is Kronecker delta function. 

To construct a directional CSWT on the sphere one must 
consider rotations on the sphere. Rotations on the sphere are 
characteristed by the elements of the rotation group SO (3), 
which we parameterise in terms of the three Euler angles 
p = (a, /?, 7). 1 The rotation of / is defined by 

[K(p)f](Lu) = fip- 1 ^, P eSO(3). 

It is also useful to characterise the rotation of a function on the 
sphere in harmonic space. The rotation of a spherical harmonic 
basis function may be represented by a sum of weighted 
harmonics of the same £ [12], [30]: 

i 

[n( P )Y im ] H = D nm(P) Ytn(Lj) , (2) 

n=-£ 

where the Wigner functions D^ n (p) are described below. 
From (2) it is trivial to show that the harmonic coefficients of 
a rotated function are related to the coefficients of the original 
function by 



[K(p)f\t 



(3) 



! We adopt the zyz Euler convention corresponding to the rotation of a 
physical body in a fixed co-ordinate system about the z, y and z axes by 7, 
(3 and a respectively. 



Note carefully the distinction between the indices of the 
Wigner functions in (2) and (3). The Wigner functions may 
be decomposed as [12], [30] 



-ima it 



W) e" 



(4) 



where the real polar d-matrix is defined by [30] 



min(£+m,^— n) 



d e mn (P) = U ("!)* 

t=max(0,m-n) 



[(£ + m)\ (£ - m)\ (I + n)\ (£ - n)\ ] 1/2 
(t + m-t)\(l-n-t)\(t + n- m)\ t\ 



cos 



2l+m-n-2t 



sin 



n— m+2t 



and the sum over t is defined so that the arguments of the 
factorials are non-negative. The Wigner functions satisfy the 
orthogonality condition 

/ dp D^ n (p) D^, n , (p) = wr — r Sw <W 5nn> , (5) 

JS0(3) 



2£+ 1 



where dp = sin (3 da d/3 dj. Recursion formulae are available 
to compute rapidly the Wigner d-matrices in the basis of either 
complex [11], [31] or real [32], [33] spherical harmonics. 

III. Analysis 

We adopt the strategy suggested by [1] and define spherical 
wavelets and the spherical dilation operator in harmonic space. 
We extend this idea here to a directional wavelet transform on 
the sphere. The dilated wavelet Sfr(u;R) is defined in harmonic 
space by 



21 + 1 
8tt 2 



(6) 



where ^£ m (R) are the spherical harmonic coefficients of 
Sfr(u;R), T m (q) are the family of wavelet generating func- 
tions and R G is the real, strictly positive dilation 
parameter. Note that the wavelet generating functions are 
not defined on the sphere but rather on the non-negative 
real line: T m G L 2 (R + , dx), m G Z (although we consider 
only \m\ < £ for \E^ m (,R), we are in general free to define 
T m , Vm G Z). For the functions Sfr(u;R) to classify as 
wavelets they must satisfy certain admissibility criteria so 
that the analysed function may be reconstructed perfectly 
from its wavelet coefficients. We consider admissibility and 
the synthesis of a function from its wavelet coefficients in 
section IV. 

An overcomplete wavelet basis on the sphere may be 
constructed from the following spherical wavelet family: 

{[n( P m(u;-R) |peSO(3), ReRt}. 

The directional CSWT of / G L 2 (S 2 , dfi) is given by the 
projection onto each wavelet basis function in the usual 
manner: 

W$(R,p) = [ dn(co) f(cj) [Kip^n^R) . (7) 

The transform is general in the sense that all orientations 
in the rotation group SO (3) are considered, thus directional 
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structure is naturally incorporated. It is important to note, 
however, that only local directions make any sense on S 2 . 
There is no global way of defining directions on the sphere 2 - 
there will always be some singular point where the definition 
fails. 

It is useful to represent the CSWT in terms of harmonic 
coefficients. Substituting into (7) the spherical harmonic ex- 
pansion of the analysed function and of the rotated wavelet, 
and noting the orthogonality of the spherical harmonics given 
by (1), one may represent the CSWT by 

oo I £ 

wUr, p) = E E E D ™(p) n n (R) ■ (8) 

£=0 rn=-ln=-t 

IV. Synthesis 

For the multi-resolution analysis we describe in the pre- 
vious section to classify as a wavelet analysis, it must be 
possible to reconstruct perfectly the original function from 
its wavelet coefficients and the wavelet basis functions. To 
ensure that this is possible the wavelets must satisfy certain 
admissibility criteria. We derive the admissibility criteria and 
perfect reconstruction formula in this section. We then relate 
the admissibility condition for each wavelet basis function to 
an equivalent condition for the family of wavelet generating 
functions. 

Consider the expression 

/ -B- / dpW$(R,p)[K(p)L r *](<j;R), (9) 

JO K JSO(3) 

where the Ly operator in L 2 (S 2 , dQ) is defined by the action 



(10) 



on the spherical harmonic coefficients of functions 
g G L 2 (S 2 , dft) (C{ is defined below). Substituting the 
harmonic representation of the wavelet coefficients given by 
(8) and the harmonic representation of \R{p)L^](uj\ R) 
into (9), and noting the orthogonality of the Wigner functions 
defined by (5), (9) may be rewritten as 



oo t 



£=0 m=-£ T 

rf\*en(R)f. (ID 



For perfect reconstruction we require the admissibility criteria 



< Cl 



8w 



2£- 



! c POO i T) 

T E I -^l^-(i?)| 2 <oo, (12) 



to hold W. We call the admissibility constant. If admis- 
sibility is satisfied it is apparent from (11) that the original 
function may be reconstructed by 

POO _J 7~> p 

fM= / -R / d PW& R >P) 
JO H JSO(3) 

x [K{p)L^\{w-R). (13) 

2 There is no differentiable vector field of constant norm on the sphere and 
hence no global way of defining directions. 



We relate the admissibility criteria that each spherical wavelet 
must satisfy to an equivalent admissibility condition for the 
family of wavelet generating functions. Substituting (6) into 
(12) and making the change of variables q = £R, the 
admissibilty constant may be rewritten as 



< C{ 



dq 



\T m (q)\ 2 <oo. 



(14) 



Since J °° -y |Y n (g)| 2 is always non-negative it is possible 
to recast the admissibility condition on the family T m in the 
following form: J °° |Y m (<7)| 2 < oo, Vm G Z and 3m e Z 
such that / °° f |rJ(g)|V0. 

V. Properties 

In this section we describe some properties of the wavelet 
basis functions and the wavelet coefficients of the CSWT 
derived in the previous sections. 

A. Compensation 

Admissible wavelets are compensated: 



/ 

Js 2 



dQ(co)^(u;;R) =0, Vi2 G 



We show this property by noting that admissibility imposes 
the condition Y m (0) = 0, Vm G Z. Considering the spherical 
harmonic transform of the wavelet, 



^M±l T m («) = J dQM R) Yf m {w) , 

and setting t = and m = we obtain the compensation 
property. 



B. Wavelet energy 

The energy of a wavelet at a particular dilation scale is given 

by 



\*(R)f = [ dn(u)\*(cu;R)f 
Js 2 

oo I 

= E E i*u;?)i 2 



£=0 m=-£ 
oo £ 



EES 1 iT m (^)i 2 



=0 m=-£ 



8tt 2 



where in the second line we have made use of the orthogo- 
nality of the spherical harmonics described by (1). 

C. Function energy 

The energy of the analysed function is given by 



[ dOH|/H| 

Js 2 

oo £ 

EEi/< 



I 2 



£=Q \ 
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where in the second line we have made use of the orthogo- 
nality of the spherical harmonics described by (1). We may 
also relate the energy of the function directly to its spherical 
wavelet coefficients: 

POO TJD P 

\\ff= d P W^ f (R,p)W(*(R,p) . 

JO K ^SO(3) 

To show this property we substitute the harmonic represen- 
tation of the wavelet coefficients given by (8) and note the 
orthogonality of both the spherical harmonics and Wigner 
functions described by (1) and (5) respectively. 

D. Wavelet domain variance 

The variance of the wavelet coefficients for a particular scale 
is used by [1] as a measure of the concentration properties of 
the wavelet transform. However, this measures the variance 
of the wavelet coefficient values and not the localisation of 
energy in wavelet space. In any case, it is useful to relate 
variance in the wavelet domain to the energy of the analysed 
function and wavelet. 

The variance of the wavelet coefficients for a particular scale 
is given by 



where we define the expected value over rotations of 

G(p) e L 2 (SO(3), dp) by 

(G(P)) = 7^[ d P G(p). 



8tt 2 



SO(3) 



The expected value of the wavelet coefficients is given by 



wUr,p)) = Y1 E E hm*UR) (C(p)> 



£=0 1 



and 



* 2 L 

Jo 



dae 1 



f 

Jo 



d/3 sin/3d£ n (/?) 



d 7 e m7 
l 



- $m0 $n0 



= 5m ^m0 ^nO j 



Pe(x) dx 



where we have applied the Wigner D-function decomposition 
described by (4). The expected value of the wavelet coeffi- 
cients is thus 



which is zero for an admissible wavelet. The variance of the 
wavelet coefficients is therefore given by 

AW{(R)] 2 = 

oo £ £ 1 

\fe m \ 2 \^M\ 2 , 



= (jw{(R,p)\ J 

oo £ £ 

= 2£ 



£=0 m=-£i 



1 



where the second line follows from the harmonic represen- 
tation of the wavelet coefficients given by (8) and by the 
orthogonality of the Wigner functions described by (5). 



VI. Wavelet construction 

In the spherical wavelet framework outlined herein wavelets 
on the sphere are constructed directly in spherical harmonic 
space. Spherical wavelets may be constructed in harmonic 
space from the analogue of the Fourier definition of planar 
Euclidean wavelets. We must then check that the candidate 
wavelets are admissible. In this section we demonstrate the 
construction of a directional spherical wavelet with an exam- 
ple. 

We construct the spherical analogue of the Morlet wavelet 
using our spherical wavelet formalism. Note that this differs 
from the definition of the spherical Morlet wavelet constructed 
from the projection of the planar Morlet wavelet on to the 
sphere [24], [25]. The planar Morlet wavelet is defined by a 
Gaussian in Fourier space centered on the wave vector of the 
wavelet, hence we define the analogue of the Morlet wavelet 
on the sphere by a Gaussian in spherical harmonic space. 
Various associations between Fouier and spherical harmonic 
space are possible. Here we choose to assocate £ and m in 
spherical harmonic space with the x and y components of a 
vector in Fourier space respectively. 3 The wavelet generating 
functions are defined by 



r m (£R) = e 



(Ifl- L) z + (m- MY 



(£R) z +L z + (m-M) z 

e ^ 5 (15) 



where L e N and M e Z, \M\ < L define the centre of 
the Gaussian when R = 1. The correction term subtracted 
is included to ensure admissibility (and takes a similar form 
to the correction used to ensure the planar Morlet wavelet 
is admissible [25]). For numerical purposes, when L > 5 
the correction is numerically negligible and may simply be 
ignored. Moreover, for numerical purposes it is possible to 
relate the maximum spherical harmonic ^ max that must be 
considered to L and R: ^ max = \(L + \J — 2 \ne)/R], where 
e is the allowable relative error and \x] is the smallest 
integer greater than or equal to x. We adopt the relation 
4ax = \{L + 5)/if|, which ensures e < 10 -5 . 

Before proceeding it is necessary to check that the candidate 
wavelet generating functions generate admissible spherical 
wavelets. We must show that he admissibility integral for the 
wavelet generating functions, given by 



f 

Jo 



dq 



T m ( 9 )| S 



-[L 2 + (m-M) 2 



dq 



e- q "(e qL -l)' z . 



converges. To do this, we split the range of integration up 
into the interval (0, e), where e > is small, and the interval 
(e, oo). For the latter, we can appeal to the comparison test 
for integrals since the ratio of -e~ q (e qL — l) 2 to q~ 2 tends 
to zero as q — > oo, for any L, and J e °° |f is finite (provided 
e > 0). For the lower range, we note that the expansion of the 
integrand for small q is qL 2 , and J^dqqL 2 is finite. 

In this section we have demonstrated the construction of 
admissible wavelets on the sphere in the spherical wavelet 
framework described herein. In particular, we have constructed 

3 Alternatively one may choose to associate t with the length of a vector in 
Fourier space and m with its phase. We have also considered this association, 
although the support of the resulting functions is not so compact. Moreover, 
we have not yet shown that the candidate wavelets that result are admissible. 
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the analogue of the Morlet planar wavelet. In Fig. 1 and 
Fig. 2 we plot this spherical wavelet for (L,M) = (10,0) 
for a dilation scale of R = 1 and R = 0.5 respectively. 
Note that this construction yields a spherical wavelet that is 
quite different to the spherical Morlet wavelet constructed by 
projecting the planar Morlet wavelet on to the sphere (to make 
a comparison see [24], [25]). Other wavelets on the sphere may 
be constructed in our framework in a similar manner. Although 
we have constructed complex spherical wavelets, real wavelets 
may be constructed simply by ensuring the reality condition 
T_ m (g) = (-l) m r* m (q) is satisfied. 

VII. Concluding remarks 

We have presented the construction of a new directional 
continuous wavelet transform on the sphere. Dilation on the 
sphere is performed by a scaling in harmonic space, as recently 
proposed by [1]. We extend the spherical wavelet analysis 
presented by [1] to a more general directional framework, 
allowing one to probe oriented structure of the analysed 
function. Our framework has the advantage that all wavelets 
and operators are defined directly on the sphere. Once the 
wavelet analysis is performed, the original function may be 
reconstructed perfectly from its wavelet coefficients and the 
wavelet basis functions. We have proved the corresponding 
synthesis formula and have derived the resulting admissi- 
bility criteria that the wavelets must satisfy. Wavelets are 
constructed in our framework in spherical harmonic space. We 
have demonstrated the construction of a wavelet basis with 
the example of the analogue of the Morlet wavelet in our 
spherical wavelet setting. Other wavelets on the sphere may 
be constructed in our framework in a similar manner. 

Not only is the extension of a wavelet analysis to a spher- 
ical manifold of important theoretical interest, it is also of 
important practical use in order to extend the potential benefits 
afforded by a wavelet analysis to data that are defined on the 
sphere. For data sets of practical size, that may contain tens of 
millions of pixels, it is imperative to have a fast algorithm for 
performing the spherical wavelet analysis. In a recent work 
we have presented fast algorithms to perform a directional 
CSWT [34] (for the CSWT derived by [22]-[28]). Fortunately, 
for a given dilation the analysis formula described herein is 
identical to that considered in these works. Our fast algorithms 
are therefore directly applicable to our new directional wavelet 
construction. 

The application of our spherical wavelet transform to prac- 
tical data is beyond the scope of this paper. In a future work 
we intend to examine the localisation properties of spherical 
wavelets constructed in our framework using numerical sim- 
ulations and to apply these to analyse the anisotropies of the 
cosmic microwave background (CMB). Spherical wavelets are 
a powerful tool for analysing the CMB, from which we are 
able to learn a great deal about the physics of the early universe 
[35]-[37]. 
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